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Abstract 



Extended Abelian monopoles are investigated in SU(2) lattice gauge theory in 
three dimensions. Monopoles are computed by Abelian projection in several 
gauges, including the maximal Abelian gauge. The number N m of extended 
monopoles in a cube of size m 3 (in lattice units) is defined as the number of ele- 
mentary (l 3 ) monopoles minus antimonopoles in the cube (m = 1,2,.. .). The 
distribution of l 3 monopoles in the nonlocal maximal Abelian gauge is shown 
to be essentially random, while nonscaling of the density of l 3 monopoles in 
some local gauges, which has been previously observed, is shown to be mainly 
due to strong short-distance correlations. The density of extended monopoles 
in local gauges is studied as a function of (3 for monopoles of fixed physical 
"size" (m/P = fixed); the degree of scale violation is found to decrease sub- 
stantially as the monopole size is increased. The possibility therefore remains 
that long distance properties of monopoles in local gauges may be relevant to 
continuum physics, such as confinement. 
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I. INTRODUCTION 



The dual Meissner effect was suggested as possible mechanism for confinement in non- 
Abelian gauge theories more than a decade ago by Mandelstam |J and t' Hooft @||, and has 
recently been the subject of intensive lattice investigations PHT7I. In this picture of confine- 



ment, (color) magnetic monopoles are presumed to condense in the vacuum, forcing (color) 
electric fields between two sources to be squeezed into a flux tube. Magnetic monopoles 
are well understood to result in confinement in compact Quantum Electrodynamics, as has 
been demonstrated both analytically f|T8 , 19f1 and numerically |20|j2T[] . The conjecture is that 



such "Abelian" degrees of freedom also drive confinement in nonAbelian theories, 't Hooft 
proposed that this mechanism can be realized in nonAbelian theories by imposing a gauge 
fixing condition that is invariant under a Cartan subgroup U^l)^ -1 of the original SU(iV) 
theory, referred to as an Abelian projection ||. Singularities in the gauge-fixing condition 
are identified with monopole world lines in four dimensions (d = 4) and monopole points 
(instantons) in three dimensions (d = 3). 

A lattice implementation of Abelian projection was formulated in Ref. in which 

several gauge- fixing conditions were also developed (following 't Hooft 0) that have been 
widely studied. The so-called maximal Abelian gauge is defined by gauge transformations G 
which maximize the following quantity |JJ (we henceforth restrict our attention to SU(2)): 

R = J2 Tr(a 3 G(x)U(x, {x)G\x + p)a z G{x)U\x, fi)G ] {x + fi)) , (1) 

where U(x, (x) are link variables. In the continuum limit, max(i?) reduces to the renormal- 
izable differential gauge (<9 M ± igA^A^ = 0, where A^ = (A 1 ^ ± iA 2 ^)/y/2. 

Local (generally nonrenormalizable) gauges can be defined by the diagonalization of an 
adjoint operator $ ||: 

- G{x)®{x)G\x) = ( ^ J a(x) ) . (2) 

Examples of Eq. (0) are diagonalization of a plaquette or a Polyakov line 0] . We will refer 
to a gauge as "local" when the gauge condition can be imposed on a site-by-site basis. 

Notice that Eqs. ([I]) and (|2|) are invariant under a local U(l) transformation G — > dG, 
where d is a diagonal matrix d = cos 9 + ia^ sin 9 |22|]. The Abelian projection of a link U is 



defined, after gauge-fixing, as its component u in the corresponding subspace [|,|5[] 

G(x)U(x, jx)G\x + fi) = w(x, fi) u(x, fl), Tr(a 3 w) = (3) 

(w^w = v)u = 1). The Abelian projection u(x,fi) can be written as u(x,fi) = 
exp[icr 3 0(x, fi)]. Monopoles are defined from the phases (as described in Sec. II) and 
are associated with singularities in the gauge fixing which, in the context of Eq. (0), occur 
at points where a(x) =0 0. 

One of the main issues that has been addressed in lattice investigations is whether the 
Abelian monopole density in a particular gauge exhibits scaling (thus having a well defined 
continuum limit). SU(2) lattice simulations in maximal Abelian gauge show scaling behavior 
for monopoles defined on elementary cubes of the lattice in d — 3 |l~2]l, and there is some 
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evidence for scaling in d = 4 PJlQl. On the other hand, it has been well established that the 
density of elementary monopoles in a variety of local gauges does not scale (this includes 
diagonalization of a plaquette and the Polyakov loop) fT0|j . 

As a result, the maximal Abelian gauge has become widely regarded as the only known 
gauge which may yet establish the role (if any) of Abelian monopoles in confinement. How- 
ever, there is as yet no compelling reason why the degrees of freedom responsible for con- 
finement in a gauge-invariant theory should only be manifest in this particular gauge. In 
fact, some monopole operators have recently been constructed which show evidence in lo- 
cal gauges for a spontaneous symmetry breaking that is correlated with the deconfinement 
phase transition fT6| , [TT|l . 

We reconsider the lattice scale dependence of the monopole density in local gauges by 
analyzing the properties of extended monopoles. We work in three-dimensional SU(2) lattice 
gauge theory. The number N m of extended monopoles in a cube of size m 3 (in lattice units) 
is defined as the number of elementary (l 3 ) monopoles minus antimonopoles in the cube 
(m = 1,2,...). A physical reason for studying extended monopoles is that confinement 
occurs on some finite physical length scale. The nonscaling behavior that has been observed 
in some gauges for elementary monopoles might therefore not be a good criterion for ruling 
out contact with the physics relevant to confinement. 

Indeed, we find that nonscaling of the monopole density in some local gauges is due 
mainly to short-distance monopole-antimonopole fluctuations, which diverge in number as 
the coupling ft — > oo. On the other hand, such fluctuations are essentially absent in the 
maximal Abelian gauge. This is reflected in our calculation of the average minimum sepa- 
ration between a monopole and the nearest neighboring antimonopole, (r m ; n ). In maximal 
Abelian gauge (r min ) scales with ft, while it vanishes (in physical units) with increasing ft 
in the local gauges that we considered. [The dimensionless quantity ft = 4/(g 2 a), where 
the coupling constant g has units of (mass) 1 / 2 in three dimensions. We work throughout in 
lattice units, where the spacing a = 1.] 

Correlations between elementary monopoles and antimonopoles as functions of their sep- 
aration are also calculated. Strong short-distance correlations are found in local gauges. On 
the other hand, in maximal Abelian gauge the monopole distribution is essentially random 
(the monopoles form a "plasma"). We speculate that large short-distance fluctuations in the 
gauge fields lead to the divergence in the number of monopoles in local gauges, while in the 
maximal Abelian gauge, which is imposed in nonlocal way, these fluctuations are effectively 
smoothed out. 

We then consider the density of extended monopoles. In local gauges the density is 
found to decrease rapidly as the monopole "size" m is increased, the rate of change being 
maximal for m near unity. This can again be understood as due to an averaging over 
short distance fluctuations (i.e., as an averaging of the charges of elementary monopole- 
antimonopole pairs over larger volumes). By contrast, the density varies slowly with m 
in maximal Abelian gauge. We also consider the density of extended monopoles of fixed 
physical "size" (m/ft = fixed), as a function of ft. For a finite value olm/ft the density of 
monopoles in a typical local gauge does not scale; however, we find that the degree of scaling 
violation decreases substantially as the physical monopole "size" increases. 

These results support the argument that strong short- distance fluctuations in the gauge 
fields can be effectively averaged out by a nonlocal gauge condition, as in the maximal 
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Abelian gauge, or (to a large extent) by considering extended monopoles on sufficiently long 
length scales in local gauges. 

The organization of the rest of this paper is as follows. In Sec. II the various gauge fixing 
schemes and observables considered in this work are described. Our results are presented in 
Sec. Ill, and in Sec. IV we give our conclusions. 



II. METHOD 

Abelian projection was carried out in the maximal Abelian gauge and in a number of local 
gauges. The maximization of R, Eq. (|l|), which defines the maximal Abelian gauge must 
be implemented iteratively. Following Ref. || we repeatedly sweep through the lattice, 
maximizing R locally by solving for G(x) analytically at each site, keeping G(x + /a) at 
neighboring sites fixed. These iterations are repeated until G at all sites becomes sufficiently 
close to the identity: 

max{l - |Tr G(x)} < 5 < 1, (4) 

with 5 = O(10~ 7 ) used as a stopping criterion. Note that this iteration procedure is not 
guaranteed to reach the the global maximum of R ||. Moreover, the iteration procedure 
is not guaranteed to increase TrG(x) monotonically over successive iterations; hence the 
stopping criterion Eq. (|]) is ambiguous. As a result, in practice the monopole number is 
not uniquely defined in the maximal Abelian gauge. For example, two field configurations 
differing only by an SU(2) gauge transformation will not, in practice, necessarily lead to the 
same number of Abelian monopoles in an actual calculation in maximal Abelian gauge. The 
implementation of the maximal Abelian gauge fixing also dominates the computational cost 
of the simulation, due to the large number of iterations (typically 1000) required to satisfy 
Eq. (D. 

A local gauge that we consider is a simple implementation of Eq. (^): diagonalization 
of the field strength tensor, $ = F12 (cf. Ref. [|J). Note that Abelian projection in this 
gauge has the property that a pure-gauge configuration (a gauge transformation of links 
equal to the identity) will in general lead to a nonzero monopole density. Furthermore, the 
monopole numbers extracted from gauge-equivalent configurations in this gauge will not be 
equal in general. Calculations were also carried out in gauges defined by the diagonalization 
of rectangular Wilson loops of various sizes, as well as in an "unfixed" gauge (defined by 
$ = 1 in Eq. (H)). These calculations lead to the same conclusions as drawn from the results 
presented here for $ = F12. 

Following Ref. |4[] we compute phases of the Abelian projections of the gauge-fixed links 
u(x, fx) = u (x, jj) + ia • u(x, /i) (cf. Eq. ©): 



<J)(x, fi) = tan 1 



u 3 (x,ft) 



_u (x,£l) 

Reduced plaquette angles <fi are then defined according to pOfl : 



x,fl) e (-7T,+7r]. (5) 



2tt7V s , 0G(-7r,+7r], (6) 
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where N s is identified with the number of Dirac strings passing through the plaquette (N s G 
[—2,2]). The number of "elementary" monopoles N m= i contained in a cube of size l 3 is 
equal to the sum of Dirac strings N s passing through the oriented lxl plaquettes on the 
surfaces of the cube. 

We define the number of extended monopoles N m in a cube of size m 3 as the number 
of elementary monopoles minus antimonopoles in the cube. The density p m of extended 
monopoles of "size" m is given by 

Pm = ^ E \N m (xi)\, (7) 

where Xi labels the coordinates of the cubes of size m 3 . We take m to divide the lattice 
length L (although this requirement can be relaxed); hence the total number of extended 
monopoles on the lattice vanishes. 

The definition of N m above corresponds to "type-II" extended monopoles defined in Ref. 
@; a "type-I" extended monopole number was also defined, which is computed from the 
phases ofmxm plaquettes on the surface of a cube of size m 3 (using a suitable generalization 
of Eq. (|(|)). In an earlier version of this work we considered type-I extended monopoles. 
However, for large m the number of type-I monopoles approaches the strong coupling limit, 
since the phases of the m x m plaquettes become essentially random [^^J. Since we are 
interested in weak coupling continuum physics we consider only type-II monopoles in the 
rest of this paper. 

The simulations were done mainly on a 24 3 lattice, using a bath algorithm at several 
values of (3 in the scaling region. The string tension on lattices of this size is found to scale 
for p > 6 (see, e.g., Ref. [p5[|). Some data was also taken at smaller values of (3 for compar- 
ison. We find that finite volume effects on the 24 3 lattice, as measured by the expectation 
value of the Polyakov line, become noticeable for (3 ^ 12 (which is the largest value in our 
data). Measurements were made on an ensemble of 500 configurations each separated by 
100 updates, which is significantly longer than the autocorrelation time for any observable 
that we considered (the maximal Abelian gauge exhibits the longest autocorrelation time of 
the various gauges). We also used 200 configurations on a 36 3 lattice to obtain some results 
at large m (these data are identified explicitly in the following). 

III. RESULTS 

The density p\ of elementary (m = 1) monopoles is shown in Fig. 1 in maximal Abelian 
gauge (cf. Eq. ([I])) and Fi 2 -gauge ($ = F i2 in Eq. (§)). It is clear that only maximal Abelian 
gauge exhibits scaling in pi, and that the difference between the two gauges becomes more 
pronounced as j3 is increased. 

Notice that maximal Abelian gauge has far fewer monopoles than the local gauges. 
To help understand this difference it is instructive to look at the spatial distribution of 
monopoles and antimonopoles. This is done by counting the number of monopoles N(r min ) 
for which the nearest antimonopole is a distance r min away. The results are shown in Fig. 2 for 
the Fi 2 -gauge. As expected with a relatively large number of monopoles and antimonopoles, 
there are many short distance pairs; in fact, most monopoles have an antimonopole in a 
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neighboring lattice cell. A natural question is whether the monopoles and antimonopoles 
are correlated. Some insight into this question is obtained by comparing the result of the 
lattice simulation with iV(r min ) for a completely random distribution of monopoles and 
antimonopoles. This is also shown in Fig. 2. A random distribution leads to an N(r min ) 
which is much broader than that observed in Fi 2 -gauge. This is indicative of a strong short- 
range correlation between the monopoles and antimonopoles in this gauge, with significant 
enhancement of nearest-neighbor pairs. 

Figure 3 shows iV(r min ) for the maximal Abelian gauge at (3 — 8. Here the number of 
monopoles with nearby antimonopoles is very small so statistical errors are large, but iV(r min ) 
is roughly consistent with a random distribution. There certainly is no enhancement of pairs 
with small separations as in Fi 2 -gauge. 

Figure 4 shows the average value of the minimum separation of a monopole and anti- 
monopole (r min ) for the two gauges as function of the number of monopoles (corresponding 
to (3 in the range 5 to 12). Again, for comparison, the results expected for a completely 
random distribution are plotted. For monopoles in maximal Abelian gauge we find rea- 
sonable consistency with a random distribution (a monopole "plasma") but in Fi 2 -gauge 
significantly smaller values of (r m i n ) are observed. 

If the monopole density scales and the distribution is random we certainly expect (r min ) 
to scale with (3. This is shown in Fig. 5. For maximal Abelian gauge, (r min ) scales very well. 
In F 12 -gauge (r min ) vanishes in physical units as (3 increases. 

These results support the idea that the excess of monopoles found in local gauges relative 
to maximal Abelian gauge is due to highly-correlated short distance fluctuations. This is 
further illustrated by considering the density p m of extended monopoles. 

In Fig. 6 we compare the density p m in the two gauges at (3 — 6. The density in 
maximal Abelian gauge decreases slowly with to; the decrease becomes more pronounced 
for to around the value of (r min ) for l 3 monopoles in that gauge. This is consistent with the 
fact that strong short-distance correlations are absent in maximal Abelian gauge, with l 3 
monopoles forming a plasma, distributed on a length scale of 0((r m i n )). By contrast, the 
density in F 12 -gauge decreases rapidly with to, with the rate of change being maximal near 
to — 1. 

Figure 7(a) shows that the density of extended monopoles of fixed physical "size" (to/ (3 = 
fixed) scales well in maximal Abelian gauge. In Fi 2 -gauge the density does not scale for any 
finite value of m//3, as illustrated in Fig. 7(b). The density p(3 3 in physical units increases 
roughly linearly with (3 for fixed m/(3. However, the degree of scale violation decreases 
substantially as the physical monopole "size" increases. This is made evident by a plot of 
the slope A(p(3 3 ) / A/3 as a function of m/(3, given in Fig. 8. This is again consistent with a 
cancellation of short- distance monopole-antimonopole fluctuations. 

IV. CONCLUSIONS 

The results presented here demonstrate that the nonscaling of the density of elementary 
(l 3 ) monopoles in a variety of local gauges, which has been previously observed, is due 
in large part to strong short distance correlations in those gauges. The distribution of l 3 
monopoles in maximal Abelian gauge was shown to be essentially random. The degree of 
scale violation in the density of extended monopoles in local gauges was found to decrease 
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substantially as the monopole "size" is increased. Our results support the argument that 
strong short-distance fluctuations in the gauge fields can be effectively averaged out by 
a nonlocal gauge condition, as in the maximal Abelian gauge, or (to a large extent) by 
considering extended monopoles on sufficiently long physical length scales in local gauges. 
The possibility therefore remains that long distance properties of monopoles in local gauges 
may be relevant to continuum physics, such as confinement. 
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FIGURES 



FIG. 1. Density of elementary (m = 1) monopoles pi/3 3 in physical units as a function of (3 in 
the maximal Abelian gauge (full circles) and Fi2-gauge (full triangles). 

FIG. 2. The number of elementary monopoles iV(r m i n ) versus r m ; n at (3 = 8 in i ? i2-gauge (full 
triangles). Also shown is N(r m i n ) for a completely random distribution (open circles). 

FIG. 3. The number of elementary monopoles N( r m j n ) versus ^min £it (3 — 8 in maximal Abelian 
gauge (full circles). Also shown is N(r m \ n ) for a completely random distribution (open circles). 

FIG. 4. The average minimum monopole-antimonopole separation (r m [ n ) versus the number of 
m = 1 monopoles iVi in maximal Abelian gauge (full circles) and i<i2-gauge (full triangles). Also 
shown is (r m j n ) for a completely random distribution (solid line). 

FIG. 5. The average minimum monopole-antimonopole separation (r m i n ) as a function of (3, 
for m = 1 monopoles in the maximal Abelian gauge (main figure) and in the i ? i2-gauge (inset). 
The solid line in the inset shows a fit (r m j n )//3 = c/(3 to the data in Fi2-gauge. 

FIG. 6. Extended monopole density p m at (3 = 6 as a function of the "size" m in the maximal 
Abelian gauge (full circles) and Fi2-gauge (full triangles). 

FIG. 7. Extended monopole density p as a function of (3, for several fixed monopole "sizes" in 
physical units (raj (3 = fixed), in (a) maximal Abelian gauge and (b) Fi2-gauge. The results for 
F\2 gauge are taken from a 36 3 lattice. Straight lines are shown in Fig. 7(b) to guide the eye. 

FIG. 8. Slope A(pP 3 ) / A(3 of the extended monopole density in Fi2-gauge, as a function of 
m/(3. The slope is estimated from the data in Fig. 7(b) using the two highest available (3 values. 
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